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ABSTRACT 

High resolution, deep imaging surveys are instrumental in setting constraints on 
semi-analytical structure formation models in Cold Dark Matter (CDM) cosmologies. 
We show here that the lack of unresolved B-band "dropouts" with V > 25 mag in 
the Hubble Deep Field (HDF) appears to be inconsistent with the expected number 
of quasars if massive black holes form with a constant universal efficiency in all CDM 
halos. To reconcile the models with the data, a mechanism is needed that suppresses 
the formation of qucisars in halos with circular velocities Vqi^c 

^ 50 — 75 km s ^ . This 
feedback naturally arises due to the photoionization heating of the gas by the UV 
background. We consider several alternative effects that would help reduce the quasar 
number counts, and find that these can not alone account for the observed lack of 
detections. If reddening by dust can be neglected at early epochs, consistency with the 
optical data also requires that the luminous extent of dwarf galaxies at high redshifts 
be larger than a few percent of the virial radii of their dark matter halos, in order not 
to overpredict the number of point-like B-band dropouts. Future deep observations 
in the J and H bands with NICMOS might reveal several objects per field or 

provide even stronger constraints on the models than existing B,V, and I data. 
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1. Introduction 

One of the primary goals of models for structure formation in the Universe is to explore 
the formation process of luminous objects such as quasars or galaxies, and to explain the shape 
and evolution of their observed luminosity functions (LFs). Recent theoretical work has focused 
on the formation of galaxies within the framework of hierarchical structure formation in CDM 
cosmologies, employing either three-dimensional numerical simulations (Gnedin & Ostriker 
1997; Zheng et al. 1997), or semi-analytic modeling (Kauffman, Guiderdoni &: White 1994; 
Cole et al. 1994; Baugh et al. 1997; Somerville & Primack 1998; Haiman & Loeb 1997, hereafter 
HL97). These calculations have been quite successful in reproducing the observed Tully-Fisher 
relation in the local Universe, the faint galaxy LFs out to moderately high redshifts z ^ 3, and the 
recently discovered clustering of Lyman-break galaxies at z ^ 3 (Baugh et al. 1997; Governato et 
al. 1998; Stcidcl et al. 1998). Hierarchical models were also applied to the more complex problem 
of quasar formation (Efstathiou &: Rccs 1988; Hachnclt & Rees 1993; Small & Blandford 1992), 
and were shown to reproduce the observed evolution in the quasar LF at redshifts z 2 under a 
minimal set of simplifying assumptions (Haiman &; Loeb 1998b, hereafter HL98; see also Haehnelt, 
Natarajan & Rees 1998). 

The above models were not tested in detail at redshifts much higher than z ~ 3 because of 
the lack of related observational data. Such tests are desirable from the theoretical point of view, 
since the model predictions arc expected to be more robust at higher redshifts when the stars and 
quasars first condensed out of the well-defined set of initial conditions in the early Universe. The 
exponentially low abundance of collapsed objects at very early times reduces the probability for 
complicated merger events. Moreover, the chemical and thermal conditions inside objects which 
assembled out of the pristine intergalactic gas are better defined before metal enrichment and 
thermal feedback from multiple generations of stars started to affect the appearance of quasars 
and galaxies. To date, the deepest optical image of the Universe is the Hubble Deep Field (HDF, 
Williams et al. 1996), which surveyed a 2.3' x 2.3' field in four bandpasses down to 29 magnitude. 
Although this image has revealed a large number of galaxies of different morphological types, a 
surprisingly large fraction of the sky was found to remain dark at the HDF detection limit (e.g., 
Vogeley 1997). The Next Generation Space Telescope (NGST), scheduled for launch in 2007, is 
planned to have a sensitivity of 1 nJy in the 1-3.5/xm range, and should be able to detect globular 
cluster scale objects out to 2; ~ 10 (see Smith & Koratkar 1998). 

The properties of faint extended sources in the HDF (Madau et al. 1996) have been well 
described by detailed semi-analytic galaxy formation models (Baugh et al. 1997). On the other 
hand, the HDF has revealed only a handful of faint unresolved sources (Mendez et al. 1996), 
and none with the colors expected for high redshift objects. The purpose of the present paper 
is to show that the lack of point-like sources in the HDF can be used to constrain models for 
the formation of cosmic structure at redshifts higher than previously possible {z ^ 4). The 
observational data and color-selection criteria are summarized in § 2 and § 3, and the simplest 
versions of the semi-analytic CDM models are outlined in § 4. In § 5, we compute the optical 
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and infrared quasar number counts in these models, and show that the predictions from standard 
scenarios without a UV feedback appear to be inconsistent with the observed lack of high-z 
sources. We also show that models with a UV feedback are marginally consistent with the data. 
In § 6, we describe possible modifications of these models that could reduce the predicted counts. 
In § 7, we discuss the expected galaxy number counts. In § 8, we present our conclusions and 
argue that future infrared observations in the J and H bands with NICMOS may in principle be 
used to place even stronger constraints on this type of scenarios. 

2. High Redshift Quasar Candidates in the HDF 

Several surveys of unresolved objects in the HDF have appeared in the literature (Elston, 
Santiago, & Gilmore 1996; Flynn, Gould, & Bahcah 1996; Mendez et al. 1996). Mendez et al. 
(1996) were able to identify 14 point-like objects, with only 6 of them at magnitudes fainter than 

V = 26.50. The faintest sources in their sample tend to have rather blue B — V and V — I colors 
(see also Elston et al. 1996). These authors point out that although the star-galaxy separation 
software fails above V ~ 27.5, misclassification at fainter magnitudes should turn more galaxies 
into stars than vice versa, as the overwhelming majority of objects in the HDF are galaxies and 
diffuse and faint features disappear first. Hence the number of point-like sources detected down to 

V = 29 should be considered as an absolute upper limit to the number of stars or quasars observed 
in this magnitude interval. In analogy with star-forming galaxies, high redshifts (z ^ 3.5) quasars 
are expected to have very red B — V and reddish V — I colors because of continuum absorption 
and line blanketing from neutral hydrogen along the line of sight. As pointed out by Elston et al. 
(1996), however, none of the stellar objects in the HDF has the colors expected for high-z quasars. 
A recent detailed search for faint quasars in the HDF by Gonti et al. (1998) confirms the lack of a 
population of very red compact sources; no z > 3.5 quasar candidates were found down to a 50% 
completeness limit at ^ = 29. 

3. Identifying Sources at z ^ 3.5 

Before comparing our models directly to the observations, we show in Figure [l| the location of 
the nine faint ^ > 25 sources detected by Mendez et al. (1996) on a — ^ vs. V — I color-color 
diagram (empty squares). The solid lines in this figure show the predicted trajectories of model 
galaxies born at redshifts 3 < z < 4.5 and passively evolving for 10^ years after an instantaneous 
burst of star formation. The filled dots show the expected colors of quasars at the same redshifts. 
Unlike galaxies - which change their colors due to the evolution of their stellar population - we 
have assumed that the aging of quasars only changes the overall normalization, but not the shape 



^Magnitudes throughout this paper are in the AB system, i.e. mx.AB ~ —2.5 log[fx/erg s ^ cm ^ Hz ^] — 48.6. 
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of their intrinsic spectra (see §4 below for details of our model). While all z > 3.5 objects are 
expected to fall within the trapezoidal region shown in the figure, none of the Mendez et al. (1996) 
sources appears to have the colors predicted for a QSO or a galaxy at high-z. One caveat to the 
color-selection method is that some QSOs could have unusual colors. Indeed, the three sources 
that have {B — V) > 1 in the Mendez sample are found close to the border of the QSO selection 
region. However, all of these objects are rather bright {V=25, 25.5, 26.3), and if they are indeed 
QSOs, they are likely to have low {z < 3) redshifts. We note that if these three objects turned 
out to be QSOs at z > 3.5, then our models would underpredict the number of bright {V ~ 26) 
quasars, in addition to overpredicting the faint counts (see Figure ^ and discussion below). 

We note that it should be possible to identify objects at z > 4.5 in much the same way 
as the z ~ 4 galaxies are identified by the B dropout technique. Since at higher redshifts line 
blanketing and Lyman-continuum absorption are even more pronounced due to the larger number 
of intervening systems along the line of sight, a straightforward extension of the dropout technique, 
using /, J, and H data, should work fairly well. In analogy with the z ^ 3.5 B dropout galaxies, 
we expect V, I, and J dropout sources to flag the redshifts z ^ 5, 7, and 9. The B — V vs. V — I 
diagram used to identify z ~ 4 star-forming objects (Madau et al. 1996) could then be replaced by 
a I — J vs. J — H color-color diagram as a probe of the redshift range 4.5 < z < 8, as illustrated 
in Figure Unreddened galaxies and quasars at z > 4.5 are expected to be confined to the 
trapezoidal region indicated in this figure, with negligible contamination from foreground stars. 

4. Model Description 

In modeling the abundance of high-redshift sources, we assume that the formation of dark 
matter halos follows the Press-Schechter (1974) theory. Hence, the net rate of change in the 
comoving number density of halos with mass Mhaio is given by the derivative -^{clnc/dM), where 
dric/dM is the Press-Schechter mass function (comoving number density per unit mass). This rate 
includes a negative contribution from the disappearance of small halos in merger events; the net 
■^{dric/dM) becomes negative when mergers dominate, i.e for masses below a characteristic mass 
scale M^{z) at any given redshift z. Quasars are short-lived and we assume that their formation 
rate is proportional to the formation rate of new halos. The negative contribution to this rate from 
merger events should therefore be ignored. For relatively bright quasars at high redshifts (and for 
all cases in our standard ACDM model), the Press-Schechter formation rate is adequate since 
the typical halo mass Mhaio is above M^:{z). Otherwise, we make the conservative assumption 
that no quasars form in halos with a mass Mhaio < M^{z) and therefore set -^{dnd dM) = (see 
discussion below). 

After each halo collapses and virializes, its gas might condense into stars or a quasar black hole, 
provided it can cool efficiently from its initial virial temperature. Efficient cooling via collisional 
excitation of atomic hydrogen requires the halo mass to be at least Mmin ~ 1O^M0[(1 + 2;)/ll]^'^/^ 
(HL97). In halos that exceed this mass, a fraction of the gas might fragment into stars - 
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forming a dwarf galaxy with a stellar mass of -Mstar = e*-^haio- A different fraction eg of the 

gas might accrete onto a massive central black hole with a mass Mbh = eg halo and exhibit a 
low-luminosity quasar activity (Eisenstein & Loeb 1995; Loeb 1998). The efficiencies and eg 
depend on several parameters, such as the redshift, halo mass, or the initial angular momentum 
of the gas. However, the simplest assumption to make is that the efHciencies are the same for all 
halos. 

An important effect we need to consider is an external feedback from photoionization, which 
exerts its influence by heating the gas before it is able to cool and condense inside the dark matter 
potential wells. Such a feedback is an inevitable consequence of the background UV flux that 
builds up after the reionization epoch, when the cosmological HII regions have overlapped and 
the universe is transparent to the ionizing flux from each individual source. Several authors have 
discussed the consequences of this feedback (Efstathiou 1992; Quinn, Katz & Efstathiou 1996; 
Thoul &; Weinberg 1996; Navarro &; Steinmetz 1997). According to the spherically-symmetric 
simulations of Thoul k, Weinberg (1996), a UV background at the level expected from the 
proximity effect (J21 =0.1 — 1; Bajtlik, Duncan k. Ostriker 1988) completely suppresses the 
cooling and collapse of gas inside halos with circular velocities fcirc ^ 35 km s~^, and substantially 
reduces (by ~ 50%) the mass of cooled baryons in systems with v^ic ^ 50 km s~^. While this 
type of feedback has been mainly invoked to suppress the formation of dwarf galaxies (Ikeuchi 
1986; Efstathiou 1992; Babul & Rees 1992), it is also natural to expect that a similar mechanism 
may inhibit the formation, or the fueling, of quasar black holes in shallow potential wells. In our 
standard model we therefore impose a minimum circular velocity Vdrc = 50 km s~^ for the halos of 
both luminous quasars and galaxies that form after the reionization epoch. This constraint results 
in a minimum halo mass of Minin ~ lOi°M0[(l + z)/ll]-3/2, two Order of magnitudes higher than 
the cutoff mass obtained from the cooling argument. 

The value of can be derived from the requirement that the resulting stellar population 
reproduces the mctallicity observed at z ~ 3 in the Lya forest, 10"^ ^ -^Lya ^ 1O~^Z0 
(Tytler et al. 1995; Songaila & Cowie 1996; Songaila 1997). Assuming a Scalo (1986) stellar initial 
mass function (IMF), and using carbon as a measure of metallicity, this yields 0.017 ;^ ^ 0.17 
(HL97). Similarly, the value of eg for high-rcdshift quasars can be found by matching the 
observed LF of quasars at intermediate redshifts, (f){L,z). Assuming that all quasars shine with a 
universal light-curve in Eddington units [i.e. that L{t) scales linearly with the black-hole mass, 
L{t) = M\)-^f(t) = egMhaio/(i)], the constant eg and the function f{t) can be adjusted so as to 
produce the observed (j){L, z) at 2.6 ^ z ;^ 4 based on the relation 

(l>{L,z)=rr dM^,i,dz'-^p^S[L-eQM^,U{h,,^)], (1) 

Jz JMn^in aMhi,ioaz' 

where i^.z' is the time elapsed between the redshifts z and z'. To account for the decline of 
the observed quasar population at redshifts below z ~ 2.6 would require additional modeling, 
but should not effect our conclusions here. The fit to the observational data for <p{L, z) (in the 
rest-frame B band) given by Pei (1995) yields eg ~ 6 x 10~^ and a light-curve shape that is well 
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approximated by an exponential, f{t) ~ exp(— t/tg) with to = 10^-^ yr in a flat A-CDM cosmology 
with Q\ = 0.65 and a Hubble constant h = 0.65 (HL98). Although the best fit values for the 
to and eg vary with cosmology (see discussion in § 6), the adopted procedure fixes eg and to 
within a factor of two in any given cosmological model (HL98). We note that could be increased 
substantially if black holes formed in only a fraction / of all halos. For example, Eisenstein & 
Loeb (1995) have argued that only the small subset of initial seed perturbations with sufficiently 
slow initial rotation can form black holes. In order to recover the fit to the luminosity function, 
the lifetime would then need to be increased to ~ to//- 

The Press-Schechter halo formation rate depends on the cosmology and power spectrum. For 
our "standard" scenario we have adopted the concordance model of Ostrikcr &; Steinhardt (1995), 
i.e. a flat ACDM model with a slightly tilted power spectrum (QQ,O,\,O,\j,h,agii-i,n)=(0.35, 
0.65, 0.04, 0.65, 0.87, 0.96). Convenient expressions for the differential volume element, 
luminosity distance, and time-redshift relation in this model were given in terms of elliptic 
integrals by Eisenstein (1997); an accurate fltting formula for the growth function is given by 
Carroll et al. (1992). We also consider two open models without a cosmological constant, one with 
Oo=0.35 and the other with f^o = 0.2. In order to remain consistent with the power-spectrum 
normalization on both COBE and cluster scales (Viana & Liddle 1996), we change the spectral 
index to n=1.15 and 1.2 respectively in these models. 

For our quasar template spectrum, we have used the average spectrum derived from the 
observations of 47 quasars by Elvis et al. (1994). The stellar template spectra have been computed 
from the low-metallicity {Z = 4x 10~"^) population synthesis models of Bruzual Sz Chariot (1996), 
assuming an instantaneous burst of star formation and a Scalo (1986) IMF. Note that a Salpeter 
IMF produces about twice as many UV photons as the Scalo IMF, for the same level of carbon 
enrichment of the intergalactic medium (IGM). We find that adopting a Salpeter IMF would 
increase our predicted number counts down to ~ 30. In our standard model, the increase is 
negligible at F = 30 but is an order of magnitude at V = 22. 



In both the quasar and galaxy models, the number of sources located between redshifts Zmin 
and Zniax per unit solid angle, with observed visible magnitude between V and V + dV, is given by 



where dVc/dzdfl is the comoving volume element per unit redshift per unit solid angle (Eisenstein 
1997), and nc{z, V) = dN/dVdVc is the comoving number density of objects at a redshift z with 
observed magnitudes between V and V + dV. 

Based on the halo formation rate d^n/dMdz, ric is given by a sum over halos of different ages 



5. 



Number Counts 
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that exist at each redshift, 

(3) 



Mi,,Uz,z'y),z' 



where the factor dM\is,\o/dV converts the number density per unit mass interval to number density 
per unit flux interval. The V magnitude in the AB system is given by 

y = -2.5 log /y -48.60, (4) 

where the observed F-band flux z', Mhaio) from a halo at redshift z with a mass Mhaio that 
has formed at z' > z equals (in units of erg sec~^ Hz""*^ cm^^), 



fviz,Z ,Mhalo) ^ 



4TTdl,{z)'^ 



roo poo 

/ diyTvii^)ew[-T{'^,z)]j{i^il + z),t^^^,) X / dvTviv 
Jo J Uo 

(5) 

Here di^{z) is the luminosity distance; e = eg = 6 x 10~^ or e = = (ilb/^o) x (0.017 — 0.17) = 
0.002 — 0.02 is the fraction of halo mass in the form of the central black hole or stars, respectively; 
t(z^, z) is the effective optical depth associated with intervening Lya forest clouds and Lyman-limit 
systems between redshifts and z at the observed frequency v (Madau 1995); Tyiv) is the 
frequency-dependent transmission of the Vgoe filteiQ; and j{i^,tz^z') is the template luminosity per 
unit (stellar or black hole) mass for either stars or quasars, in units of erg sec^ -1 Hz-i Mq^ The 
time dependence of j(i^) is obtained from our fitted light-curve function L(t) for quasars, and from 
tabulated evolutionary tracks for stars (Schaller et al. 1992; see HL97 for details). We suppressed 
the emergent flux from all objects at emission frequencies above Lya before reionization occurs 
(assuming sudden reionization at the redshift z = 13 or 9 for stars, and z = 11.5 for quasars; see 
HL98 for details), since the Gunn-Peterson optical depth of HI at these frequencies is exceedingly 
high prior to reionization. 

The expected counts of z > 3.5 quasars in the V band are shown by the solid lines, labeled 
A, in Figure ^. The bottom panel also shows the observational data points from the B band 
(Pei 1995), converted to V counts using the shape of our template spectrum. By construction, 
the differential model number counts for quasars in the narrow redshift interval 3 < z < 3.2 
are consistent with the data between 17 < V < 22.5. The upper panel in this figure shows the 
cumulative counts of z > 3.5 quasars, and reveals that down to the V = 29 limit of the HDF 
our standard model predicts about 4 objects. Since no unresolved sources at z > 3.5 have been 
found in the HDF, a model which is consistent with the HDF data at the 2a level should predict 
no more than 2 such high-z QSOs. On the other hand, since the observations are only ~ 50% 
complete, a prediction of 4 QSOs could be allowed. Our standard quasar scenario therefore 
appears to be marginally consistent with the HDF data. Although the existing data on the LF 
only constrains the individual quasar model parameters such as eg and to to within a factor 



* Available at http: / /www. stsci.edu/ftp/instrument_news/WFPC2/Wfpc2_faq/wfpc2_phot_faq.htmj . 
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of ~ 2, the z > 3.5 number counts are predicted to a much better accuracy, since the model 
is cahbrated by simultaneously adjusting eg and to to reproduce the quasar LF at z 3. As 
a result, an overprediction of the quasar counts by a factor of 2 could still be significant. It is 
important to remember here that in our standard model, we have eliminated quasars in halos 
with Vcirc < 50 km s~^ to indicate the effect of photoionization heating. The thin dashed curves in 
Figure |3|, labeled Al, show the number counts in a model without this UV feedback, and reveal a 
prediction of ~ 8 objects to V = 29. The overprediction in this model would be rather severe, a 
factor of ~ 4. 

The expected number counts of z > 3.5 galaxies in our standard model are shown by the 
curves labeled A in Figure ^. The solid and dashed lines bracket the expected metallicity range 
of the IGM at z ~ 3, W~^Zq < Z < 1O~^Z0. Depending on the metallicity calibration, the total 
number of galaxies observable down to ^ = 29 mag is between 20 and 140. This exceeds the 
predicted quasar counts by a factor of 5 — 35, and is 10 — 100 larger than what would be allowed 
for point -sources in the HDF. However, depending on their assumed intrinsic sizes, some or most 
of the z > 3.5 galaxies could actually be resolved (see discussion below). We also note that recent 
revisions of the C/H ratio in these absorbers prefer lower metallicities, closer to 10~'^Zq (Songaila 
1997). 

6. Reducing the z > 3.5 Quasar Counts 

As the simplest scenarios with a UV feedback appear to be only marginally consistent with 
the faint V counts of quasars, the question naturally arises: how should one modify the models 
in order to reduce even further the expected number of faint point-like objects? In general, the 
theoretical number counts are sensitive to the underlying cosmology and power spectrum, the 
relation between the halo mass and the emitted luminosity, as well as the time-evolution of each 
object, including its merger history. 

The simplest explanation for the lack of faint quasars would be a more effective UV feedback 
mechanism, resulting in a higher minimum halo size for luminous objects. In Figure |3|, the lines 
labeled A2 show the number counts of quasars in a model with Vdrc > 75 km s~^. The top panel 
of this figure demonstrates that with this cutoff value of Vdrc, the z > 3.5 counts are reduced 
to 2 objects, i.e. consistent with the observed lack of detections. This required cutoff value is 
higher than the values obtained from spherical simulations (Thoul & Weinberg 1996). With a 
UV background at the level of 0.3 < J21 < 1, the cooled baryon fraction is fully suppressed 
for halos with Vdrc ^ 35 km s~^, suppressed by 50% for Vdrc ^ 50 km s~^, and unaffected for 
circular velocities above Vdrc ^ 75 km s~^. However, using 3-D simulations, Navarro &: Steinmetz 
(1997) find the UV feedback to be more effective, and their circular velocity thresholds are 
Vdvc = 80 — 200 km s~^, exceeding the required value Vdrc = 75 km s~^ derived here. 

In the absence of a photoionization feedback, another straightforward approach to reduce the 
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faint counts is to make the slope of the luminosity function shallower. If the luminosity of each 
object were steady and proportional to the halo mass, L oc Af , then the LF would have been related 
to the mass function simply by 4>{L) = dric/dM x dM/dL oc dric/dM, i.e. the LF would have had 
the same slope as the Press-Schechter mass function. At low masses dric/dM oc M"^"*"^""^^/^, 
independent of cosmology and only very weakly dependent on the power spectrum index n. By 
introducing a non-linear scaling between luminosity and halo mass, L oc M^"*"" (see discussion 
below) one could tune the LF to have a shallower slope, <f>{L) oc L^^^^/(^+"\ 

In reality, each source shines only for a limited amount of time At which is much shorter than 
the Hubble time. This results in (f>{L) At x dric/dMdt x dM/dL{t), where dM/dL oc exp(— t/to) 
for our exponential light -curve. Here, a complication arises because at masses below M^, the 
time derivative of the mass function is negative, dnd dMdt < 0. Physically, this derivative is 
the sum of a positive term (the formation rate of new halos with mass M), and a negative term 
(the disappearance rate of existing halos due to mergers), and while the formation rate exceeds 
the merger rate for M > M^, the latter dominates for M < M^,. The exact ratio between the 
formation and merger rates requires a scheme that goes beyond the Press-Schechter theory and 
describes the merging history of halos (Lacey & Cole 1993; Kauffmann &; White 1993). For the 
sake of simplicity, we will assume here that no halos form below the threshold mass and that 
no halos are destroyed by mergers above this mass, i.e. we simply set the derivative dnd dMdt 
to zero for M < M^, and to the Press-Schechter value for M > M*. Although this assumption is 
ad-hoc, it can only underestimate the number of objects. Therefore, the reduction in the number 
counts derived here should be considered as the minimum amount of change that is necessary in 
order to fit the observations. 

Under the above prescription, the LF flattens to (j){L) oc 1/L at luminosities corresponding to 
M < M*. The cosmological parameters (iloj^jHo) and the power spectrum parameters (cigh-i,™), 
flatten the number counts at ^ = 29 (with V dN / dV =const) in models where low-mass halos form 
at sufficiently high redshift. Successful models raise the threshold mass M* at 3 ^ z ^ 6 to high 
enough values so that F ~ 29 objects become associated with sub-M* halos. For example, at 
z ~ 3, the halo mass corresponding to y ~ 29 in our models is ~ 10^~^ Mq. In general, M^: is high 
in models with large values of crg^-i, n, and 0,q, and low values of A. These parameters are related 
by the normalization of the power spectrum on the COBE (X-ray cluster) scale, which requires 
cjgii-i to decrease (increase) with decreasing Qq. In Figure ^, we show two A = open models, 
in which = 0.35 and 0.2 (lines labeled Bl and B2, respectively). The values of (cgh-i and n) 
in these models which are consistent with both the COBE and cluster normalizations are (0.87, 
1.15) and (0.87, 1.2), respectively. Note that the values of to and eg also need to be readjusted in 
these cosmologies. We find {to,eQ)={W^-'^yi, lO'^-^) for = 0.35 and {to,€Q)={W^-^yr, lO'^-^) 
for = 0.2. The lines in the bottom panel labeled Bl and B2 show that the z ~ 3 number 
counts indeed flattens at 1^ ~ 23 and y ~ 24 in these models. The upper panel of this figure also 
shows that in the f^o = 0.35 model, the cumulative z > 3.5 counts at ^ = 29 are lowered to the 
required level of ~ 2 objects. While the z > 3.5 counts in the f^o = 0.2 model are also reduced. 
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the predicted number of quasars at the HDF hmit is still too large in this case (i.e. ~ 3 objects, 
making the observed lack of detections a 3a event). 

A third alternative approach to the photoionization feedback for reducing the faint z > 3.5 
counts is to make the LF decline faster at high redshifts. While the quasar LF is reasonably well 
determined at z ~ 3, it is significantly more uncertain around z ^ A (see Figs. 1 and 2 in Pei 
1995). One could exploit this uncertainty and modify the models so that the predicted number of 
z ~ 4 objects is at the low end of the allowed range. We have found that the redshift evolution 
can not be made sufficiently steep by changing the cosmology or the power spectrum, while still 
maintaining a good fit to the LF around z ~ 3. However, the evolution can be made steeper if 
the halo mass associated with a fixed luminosity is increased, since the formation rate of higher 
mass halos evolves more rapidly. To achieve this change, we introduce a simple dependence of 
luminosity on halo mass, L = eo(Mhaio/10*MQ)^+" with a > 0. A similar scaling, including an 
additional redshift dependence, L oc Mbh oc f^jj.^. oc M^^j^(l + z)^/^ has been recently considered by 
Haehnelt, Natarajan & Rees (1998). Our adopted scaling lacks the factor (1 + z)^^"^ , and therefore 
predicts a steeper redshift evolution and a slightly smaller number of z > 3.5 objects than in 
Haehnelt et al. (1998). Since most objects at F = 29 are located at z < 6, the difference in the 
number counts predicted by the two models in the relevant magnitude range is small ( ^ 15%). 

For any given value of a, the fitting procedures described in § 4 can be repeated to obtain the 
values of eo and to- Although by definition our procedure for quasars fits the LF at 2.6 < z < 4, 
an increase in the value of a yields an increase in the characteristic mass of the LF, and a decrease 
in the total number of faint z > 3.5 quasars. We have experimented with values of < a < 1, 
and found that for a ^ 0.5, the fitting procedure for quasars fails, due to the implausibly steep 
redshift evolution of the model LF already for the observed redshift interval 2.6 < z < 4. In 
other words, an acceptable fit in our models to the LF at z = 3 would result in a predicted 
z ~ 4 LF that falls below the errorbars of the observed z = 4 bin in Pei (1995). Note that much 
stronger deviations from the linear relation Mbh oc Mhaio would also be mildly inconsistent with 
observational data for local galaxies (Magorrian et al. 1997; McLeod 1997). The largest plausible 
reduction in the faint counts can therefore be achieved using a = 0.5. In this case, we obtain 
Mbh/Mhaio = 8 X lO~^(Mhaio/lO®M0)°-5 and f{t) = exp(-t/3 x lO^yr), and predict the number 
counts shown by the curves labeled C in Figure ^. The bottom panel shows that the quality of 
the fit to the z ~ 3 LF at 17 < y < 22.5 is comparable to the fit in our standard model. The top 
panel shows that relative to model A, the z > 3.5 counts decrease by a factor of 3 to ~ 3 objects. 
However, even with the largest allowed value of a, the z > 3.5 counts remain a too high (i.e. ~ 3 
objects, making the HDF a 3a event). 

A fourth effect that could contribute to the deficit of observed faint quasars is the expected 
shift in the "big blue bump" (BBB) component of the quasar spectrum to higher energies for 
the lower mass black holes which populate the Universe at high redshifts. In thin accretion disk 
models (e.g., Peterson 1997) the emission temperature, and hence the peak wavelength of the 
BBB, scales with the black hole mass as Apeak oc M^^^. In our present template spectrum, which 
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is based on data from local {z ^ 1) quasars (Elvis et al. 1994), the BBB peak is tentatively 
located at lOOOA (see HL98). Since the characteristic black hole might be smaller by several 
orders of magnitude at z ~ 4-6 than at z ^ 1, the BBB feature might shift by a factor of a few 
in wavelength. Such a shift would reduce the observed V band flux by a factor of ~ 7, and the 
rest-frame B band flux (which we use to calibrate our model) by a factor of ~ 2. In particular, if 
the BBB in high-redshift quasars is blueshifted to ^ 300 A, then the faint V quasar counts would 
be reduced by a factor of ~2, not quite sufficient to bring the number counts into agreement with 
the HDF data. 

7. Number Counts of z > 3.5 Galaxies 

The predicted number of z > 3.5 galaxies in our models at the HDF limits is a factor of 5-35 
higher than the number of quasars. However, faint galaxies might have sufficiently large angular 
sizes so as to appear extended and avoid the overprediction of point sources. 

A simple estimate shows that this might indeed be the case. To estimate what fraction of 
the z > 3.5 galaxies may have angular diameters above the 0.1" pixel size of the three wide 
field cameras (WFC), we assume that each collapsed halo has a 1/r^ density profile, with an 
overdensity of 50 at the edge of the halo relative to the average background density (equivalent 
to an average overdensity of 150), and that the luminous region of the galaxy extends out to 5% 
of the virial radius of the halo (see Navarro, Frenk & White 1995). With these assumptions, we 
find that for the low star formation efficiency (producing an average metallicity of Zi^ya = 1O~^Z0 
in the intergalactic medium), only 0.5% of the z > 3.5 objects would be unresolved at V = 29.8. 
As a result, the HDF is expected to contain no point-like galaxies at this magnitude limit. The 
unresolved fractions around the HDF sensitivity limit are significantly higher for a high star 
formation efficiency (^Lya = 1O~^Z0). In this case, a fraction equal to 2%, 9%, 18%, and 24% of 
the z > 3.5 objects would be unresolved at V = 27.5, 28.0, 28.9 and 29.8, respectively. Therefore, 
given the expected intrinsic sizes of the galaxies, the overprediction problem disappears at low 
star formation efficiencies, but is not fully resolved if the star formation efficiency is higher and 
the apparent angular size of galaxies smaller. This could be interpreted as support for a low star 
formation efficiency, and a corresponding low average metallicity of the IGM at z ~ 3. 

It is interesting to examine the sensitivity of the number counts to the underlying cosmology 
and to the relation between Mgtar and Mhaio assuming that the stars occupy a region much smaller 
than 5% of the virial radius of each halo so that most z > 3.5 galaxies are unresolved. Obviously, 
it is more difficult to reconcile the number counts of such galaxies with the HDF data than that of 
quasars. There are two reasons for this result. First, the number counts of galaxies arc higher to 
begin with. Second, the model galaxies are ~ 2 orders of magnitude less luminous than quasars, so 
that the same V magnitude corresponds to a halo mass which is larger by ~ 2 orders of magnitude. 
The differential number counts flatten at the fixed mass M*; for galaxies this mass corresponds to 
F ~ 30 at z = 4.5. 
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As an example, the curve labeled B in Figure ^ shows the galaxy counts in the Qq = 0.35 
open CDM model. The counts are somewhat less steep than in the A model (the solid curve 
labeled A). Although the flattening of the mass function at M < reduces the faint counts, the 
effect is negligible at y ~ 29. The figure also shows the results of changing the dependence of 
the star formation efficiency on the halo mass (changing a). When the star formation efficiency 
is normalized to Zi^ya = 10~^Zq, then a = 2 and eg = 3.3 x 10~^ are required for consistency 
with no detections (see the line labeled C2). Similarly, when Zi^ya = W~'^Zq is assumed, a = 2.5 
and eo = 7.3 x lO"^'^ are required (CI). As the figure shows, the slopes of the latter number 
count curves are flattened considerably. For reference, we show also the number of counts of all 
extended sources detected in the HDF in the V band. These counts have a shallow slope and 
demonstrate that at the V = 29 limit only a small fraction of all detected objects are expected to 
be at 2 > 3.5. Finally, we note that an additional reduction of the galaxy counts could also be 
caused by dust extinction. Estimates based on the Calzetti et al. (1994) extinction curves indicate 
that the number of 2 ^ 3 sources could be as much as ~ 7 times higher than currently detected 
(Pettini et al. 1997), when dust reddening is taken into account. 

Given all of the above uncertainties, we conclude that, unlike in the case of quasars, there is 
no significant discrepancy between the expected galaxy counts and the lack of faint point sources. 
This is true as long as the stars occupy a region extending out to a few percent of the virial radius 
of their host halo, and the star-formation efficiency is close to the value expected from a low 
average Lya forest metallicity, Zi^ya = W^^Zq. 

8. Conclusions 

The observed lack of faint, unresolved, and red sources in the HDF down to F ~ 29 requires 
models that are consistent with the data at the 2a level to predict no more than 2 point -like 
high-z objects to this limiting magnitude. Given that the observations are ~ 50% complete, this 
appears to be marginally consistent with the predicted number (~ 4) of high redshift quasars in 
the simplest hierarchical models, which suppress the formation of black holes in CDM halos with 
circular velocities Vdrc ^ 50 km s~^. Without a photoionization feedback, these models predict 
~ 8 quasars from z > 3.5. In models without a UV feedback, the predicted quasar counts might be 
reduced by early structure formation (so that M* ^ IO^-'^Mq hy z = 3), but only if one postulates 
that new black holes do not form in halos below the nonlinear mass-scale, M^:{z). An example 
for such a model is the = 0.35 open CDM model, with a power spectrum slope n = 1.15 and 
normalization cr^^-i = 0.87. The z > 3.5 number counts can be further reduced by introducing a 
black hole formation efficiency Mbh/Mhaio that depends on the halo mass as oc M^^. Similarly, 
a shift of the big blue bump component to shorter wavelengths in the spectra of faint quasars 
could also help reduce the predicted counts. Combined with a minimum circular velocity of 
"fcirc = 50 km s~^, either of these latter two effects could reduce the number counts to less than 
~ 2 objects. On the other hand, we find that without a photoionization feedback, neither of these 
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effects would be sufficient by itself to explain the lack of detections. 

The models predict ~ 5-35 times more z > 3.5 galaxies than quasars at y = 29. If the star 
formation rates in these galaxies are adjusted to produce an average Lya forest metallicity of 
10"22q by 2 = 3^ then more than 10% of all galaxies could appear unresolved at the 0.1" pixel 
size of the WFC, possibly leading to another discrepancy with the observed lack of detections. 
The situation changes if the star formation rate is lowered by an order of magnitude (normalized 
to produce a 1O~^Z0 average IGM metallicity at z ~ 3), since then almost all of the predicted 
z > 3.5 galaxies are expected to be extended rather than point sources. 

Future observations at longer infrared wavelengths such as the J and H infrared bands, 
would be better suited to the study of the z ^ 4.5 universe as they allow bigger volumes to be 
probed. In Figure ^, we show the expected number counts of quasars in the /, J, and H bands 
for the model which is consistent with the optical HDF data. An increasingly larger number of 
objects are predicted to be detected in the /, J, and H bands as these are affected by intervening 
absorption only at extreme redshifts. These number counts could be somewhat reduced since faint, 
compact WFC sources lose signal-to-noise when observed by NICMOS, due to the poorer angular 
resolution of the infrared images. Figure ^ shows the expected number of galaxies in the infrared 
bands. At 30 mag, the gain in observing galaxies redward of the /-band is not as significant as 
for quasars, because the intrinsic flux of unreddened galaxies, unlike those of quasars, decreases 
as a function of increasing wavelength at z ~ 4. The redshifts of either type of sources can be 
obtained using a straightforward analogy to the B dropout technique, applied to the I — J vs. 
J — H color-color diagram shown in Figure ^. Future deep observations with NICMOS in these 
bands of a field comparable in size to the HDF would be able to either reveal tens of ultra-high 
redshift quasars, or else place tighter constraints on quasar models than currently possible using 
V and / data. 
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Fig. 1. — Theoretical B — V vs. V — I color-color diagram for z > 3 objects. The solid dots 
show quasars whose intrinsic colors are assumed not to change with cosmic time, while the solid 
lines show galaxies with ages of — 10^ years. The data points [empty squares) are the nine faint 
{V > 25), unresolved objects from Mendez et al. (1996). The trapezoidal region enclosed by the 
dashed lines should selects all unreddened sources with z <:,3.5. 
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Fig. 2. — Theoretical I — J vs. J — H color-color diagram for 4.5 ^ z ^ 8 objects. The solid dots 
show the expected colors of quasars, while the solid lines show galaxies with ages — 10^ years. 
The trapezoidal region selects objects with z ^ 4.5. 
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Fig. 3. — F-band number counts for high-redshift quasars. The lower panel shows differential 
counts for 3 < z < 3.2, as well as data adapted from Pei (1995). The upper panel shows cumulative 
number counts for z > 3.5. The curves labeled A show predictions of our standard ACDM model 
with Mbh oc Mhaio- The curves labeled Al show the counts in a model with no UV feedback. The 
curves labeled A2 correspond to a model with halo circular velocities Vdrc ^ 75 km s~^. The curves 
labeled Bl and B2 correspond to open CDM models with ^0=0. 35 and 0.2, respectively, and the 
curves labeled C corresponds to a "tilted" model with cx; MM . 
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Fig. 4. — V counts for high-redshift galaxies in the simplest model with Mgtar oc Mhaio- The 
solid lines correspond to low star formation efficiencies (Mgtar = O.OlTMhaio) that produce an IGM 
metallicity of W~^Zq, and the dashed lines correspond to a star formation efficiency which is 10 
times higher. The curves labeled A correspond to our standard model. The curve labeled B shows 
the counts in the open CDM model. The curves labeled CI and C2 correspond to models with 
Mstar oc M3^i„ (for ^Lya = ^O'^^q) and M^tar oc M^fj^ (for ^Lya = IQ-^Zq), respectively. The 
solid dots show all extended sources detected in the HDF in the V band (Williams et al. 1996). 
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Fig. 5. — V, I, J, and H counts for quasars in the ACDM model with 00=0-35 and minimum 
circular velocity Vdrc = 75 km s~^. This model is consistent with HDF data in the V band at the 
5% level. Sensitivities shown for each band are for a fixed signal-to-noise ratio and exposure time. 
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Fig. 6. — Theoretical V, I, J, and H counts for z > 3.5 galaxies in our standard ACDM model 
with Zi^ya = W~^Zq. Since most galaxies at F < 30 mag may be resolved, this model could still 
be consistent with the HDF images. 



